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The present effort explores the influence of a longitudinal applied electric field on the magne¬ 
tocaloric properties of graphenes. The magnetic entropy change AS(T, AB, E) has two contribu¬ 
tions: S cos (T, B, E ) with an oscillating character on m, inversely proportional to the magnetic field; 
and S per (T, B, E) with a periodic character also on m. In comparison to the case without electric 
field, the maximum value of the magnetic entropy change either increases or decreases due to an 
applied electric field, depending on the value of m; and, in addition, the temperature in which the 
maximum entropy change occurs decreases due to the electric field. 
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The magnetocaloric effect (MCE) is an intrinsic and 
exciting property of magnetic materials, in which, under 
an isothermal process, they are able to exchange heat 
(A Q = TAS) with a thermal reservoir, under a magnetic 
field change AB. This amount of heat is measured from 
the magnetic entropy change AS(T, AB), usually used 
to characterize materials concerning the magnetocaloric 
properties. From the adiabatic process, under a magnetic 
field change, the magnetic material changes temperature. 

In addition, the MCE has a significant technological 
importance, since magnetic materials with large MCE 
values can be employed in various thermal devices 7 - 
promising candidates to substitute, in a near feature, 
the common devices, i.e., those based on compression- 
expansion of Freon-like gases (for instance, air condition¬ 
ers and household refrigerators). This probable scenario 
arises since magnetic cooling machines fit in a clean, safe 
and sustainable technology; where these systems do not 
use any noxious gases to the ozone layer and, in addi¬ 
tion, have greater efficiency of cooling power and lower 
energy consumption when compared with the traditional 
devices. However, the magnetocaloric potential of ma¬ 
terials is not limited to room temperature applications; 
and the best example is the adiabatic demagnetization 
refrigerator, that can reach rnK scale 7 . 

In spite of those direct applications, the inverse mag¬ 
netocaloric effect (negative magnetic entropy change), is 
opening doors for devices, since the reverse operation can 
be done (cooling instead of heating, when the magnetic 
material is submitted to a magnetic field change under 
adiabatic conditions). In this sense, a lot of efforts have 
been done by the community 7 7 ? ? ? ? ? 

Nowadays, the magnetic materials available and stud¬ 
ied by the scientific community do not have yet the 
needed characteristics to be used in large scale, due to 
technological and/or economical restrictions; and the lit¬ 
erature on the magnetocaloric effect sum efforts for ap¬ 
plications considering only materials with cooperative 
orderings, like ferro-, antiferro- and even ferrimagnetic 
materials 7 . However, recently, diamagnetic materials re¬ 
ceived due attention never received before 7 7 , and an os¬ 


cillatory behavior (named oscillating magnetocaloric ef¬ 
fect - OMCE), due to the crossing of the Landau levels 
through the Fermi level - , was predicted, in analogy to 
the de Haas-van Alphen effect. 

To optimize the effect, other diamagnetic materials 
have to be considered; however, 2D standard systems 
have a energy spectra similar to the 3D case and therefore 
a 2D massless Dirac systems, more precisely graphenes, 
were the natural choice 7 . 

Graphene is a planar sheet of Carbon atoms packed in 
a honeycomb lattice. Electronic properties of graphene 
are quite different from standard 3D diamagnetic mate¬ 
rials, since electrons in graphene have zero effective mass 
and these behave like relativistic particles, described by 
the Dirac equation 7 . While for a non-relativistic 3D 
diamagnetic material the Landau levels are given by a 
harmonic oscillator 7 7 : 


Ei = huj 



where 


. HeB 
nio =- 


m 


(1) 


for a 2D massless Dirac material the levels are given by 7 : 

Ej = hu:'\J~j where hid = V 2 HeBvF (2) 

where i,j = 0,1, 2, • • • represents the Landau levels and 
vp = 10 6 m/s stands for the Fermi velocity (only 300 
times smaller than the speed of light). Note the energy 
spectra for a non-relativistic material are equally spaced, 
while it is not for a 2D massless Dirac material. 

A 3D non-relativistic diamagnetic materials (like 
Gold), has the OMCE maximized at low temperatures 
(c.a. 1 K) 7 - ; while the maximum OMCE of a 2D mass¬ 
less Dirac system (a graphene), appears in a more com¬ 
fortable temperature range (c.a. 100 K) 7 , due to the 
relativistic properties of the electrons into this material, 
more precisely the huge Fermi velocity. In addition, for a 
3D non-relativistic material, the change of magnetic field 
change needed to invert the magnetocaloric effect from 
normal to inverse is c.a. 1 mT 7 7 (useful for application 
in high sensible magnetic field sensor at low tempera¬ 
tures); while for a graphene it is c.a. 3.4 T 7 . Following 
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these results, the aim of the present paper is to under¬ 
stand the influence of a longitudinal electric field on the 
OMCE of graphenes, in order to further optimize the ef¬ 
fect. 

The energy spectra for a nano-ribbon graphene sheet 
on x — y plan with E = ( —E , 0,0) and B = (0, 0, B), i.e., 
electric field along the graphene plan and magnetic field 
perpendicular to graphene plan, is ? ? ? : 

E n = fojj"y/n + hvFpky (3) 

where 

Hu>" = V2heBv F (l~ P 2 ) 3/4 (4) 


and n = 0,1, 2, • • •. Note n is positive and it means we 
are dealing with electrons, since holes (n < 0) is out of 
the scope of the present work. Above, 


P 


E 

v F B 


< 1 


( 5 ) 


and k y = 2nl/L y (l = 0,±1,±2,...) is related to the size 
L y of the graphene along the y direction. Note, increasing 
/3 the Landau levels spacing decrease; and for a critical 
value (/3 = 1), the entire Landau strucuture collapses 7 . 
From this energy spectra, the grand canonic potential 
per graphene area can be evaluated 7 7 7 7 7 7 : 


4>(T,B,B) = (6) 

MT) + 4> rao (B, E) + <f> per (T, B, E) + <f> cos (T, B , E) 

where 4>o (T) is a contribution that depends only on the 
temperature; $> no (B, E) depends on both, magnetic and 
electric fields in a non-oscillatory fashion and does not 
depend on the temperature; $ per .(T, B, E) has a periodic 
behavior and depends on the temperature, magnetic and 
electric fields and, finally, 4> cos (T, B, E) has a cosine os¬ 
cillatory character and also depends on the temperature, 
magnetic and electric fields. Thus, the entropy 

S{T,B,E) = ~^(T,B,E) (7) 

of the system can be obtained and it is easy to see that 
Sno = 0 and 

S(T, B, E) = S 0 (T ) + S per (T , B, E) + S cos {T , B, E) (8) 

As mentioned above, to characterize the magne- 
tocaloric effect, the magnetic entropy change must be 
considered 


AS(T, AB, E) = S(T, B, E) — S(T, 0, E) (9) 

However, from equation [5j we verify that the magnetic 
field can not be zero, otherwise [3 would be bigger than 
one. Thus, there is a minimum value for the magnetic 
field: B m i n = E/v f and therefore B —> B m i n implies to 
/3 —> 1. The magnetic entropy change reads then as: 

A S{T, AB, E) = [ S per (T , B, E) - S per (T, E/v F , E)] + 

( 10 ) 

[S cos {T, B, E) - S cos (T , E/v f , E)] 



FIG. 1. Periodic function A m , from equation ESI 


similar to what was find to a 3D non-relativistic 
material' and graphene without external electric field' . 

To obtain the magnetic entropy change, only two grand 
canonic potentials must be known: 4> cos (T, B, B) and 
4> per .(T, B, E). The first one, per graphene area, is 7 7 : 


<L cos (T,B,B) 


(eB) 2 v 2 F 

ny 0 


fc =1 


cos(7r km) Xk 
(-7T k) 2 sinh(xfe) 


(11 ) 

where yo = Kv F \fN o7r is the zero temperature and zero 
magnetic field chemical potential, i.e., the Fermi energy; 
and Nq = 10 16 m~ 2 is the density of charge carriers 7 7 . 
Above, k is only the index of the summation and has no 
physical meaning. In addition: 



( 12 ) 


where <j >o = 7 ih/e = 2.06 x 10 15 Tm 2 is the magnetic 
flux quantum. Finally: 


</>o k B T 1 

"b“o^ (I-/? 2 ) 3 / 4 


(13) 


where v F = hv F /2 f\/T%tt = 9.43 x 10 _38 Jm 2 . Note 
both, to and Xk are dimensionless. The other grand 
canonic potential of interest is' ' : 


%MT, B, E) = *L,Kv r ^A m Y. (M) 


where the periodic behavior lies on 

A m = [to 2 — 2m{2a + 1) + 4<t(ct + 1)] (15) 


Above, a = [to/ 2], where the [ ] brackets mean the inte¬ 
ger part of the argument. The periodic behavior of A m 
can be seen in figure |T] and it is important to note this 
function is —1 of odd values of to and, on the other hand, 
it assumes 0 for even values of m. Still considering equa¬ 
tion E3 note it depends on L x , the size of the graphene 
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along x direction. We considered L x = 10” 8 m 7 . Those 
grand canonic potentials (equations II ll and lldl) . obtained 
from references ? ? , consider all electronic levels of the 
system. 

Thus, from the grand canonic potentials above de¬ 
scribed, the corresponding entropies can be obtained: 

SUT.B.E) = 2rt B A (1 - c sAM 7{xk) 

k=l 

(16) 

and 

N °° 

S per {T,B,E) = 2TT 2 k B L x pA m — ^N^y k7(x k ) 

(17) 

Above 

7(z) = ~ \ where L(z) = coth(z) — - (18) 

smh(z) 2 

is the Langevin function. As usually considered for ther¬ 
modynamic quantities with this kind of oscillations 7 7 7 , 
we will consider only k = 1 term: the hyperbolic sine 
function at the denominator of equation[18]rapidly damp¬ 
ing the summation. Of course, it is an approximation; 
however, after this assumption, we can understand the 
physics behind this model. Thus, the relevant contribu¬ 
tions to the entropy reads as: 

S cos {T,B,E) = 2k B — (1 - /3 2 ) 3/4 cos(7rm)T(x) (19) 
m 

and 

S per (T, B, E) = 2n 2 k B L x 0A m ^-y/N^7(x) (20) 

m 

where x = X\. 

To evaluate eauationflTIl the S(T, E/vf , E ) terms must 
be considered. From equation El B ^ Bmin — E/vf 
implies to /3 —> 1 and therefore x —> oo (see equation 
m. However, it is easy to see that 7{x —> oo) —> 0; and 
therefore S(T,E/vp,E) terms also tends to zero. Thus, 
the magnetic entropy change per graphene area reads as: 

A S{T, A B, E) = S per (T, B, E) + S cos (T, B, E) (21) 

where the above contributions are in equations [T5] and 

m 

Considering first the case without electric field, i.e., 
0 = 0 (discussed in reference ? ), S per {T, B, E) 

term vanishes and therefore S cos {T , B, E) term promotes 
oscillations on m. It is represented in figure [2] as 
AS(T, AB, 0). To understand the influence of the elec¬ 
tric field on the magnetocaloric properties of graphenes, 
0 must be close to 1/2 (i.e. 2 E s=s vfB)\ otherwise, ei¬ 
ther S cos (T, B, E) » S per (T,B,E) or S cos (T,B,E) <C 
S per (T, B, E). Note, from the case without electric field, 
the oscillations rapidly damping and for m > 4 there 
are no longer significant oscillations. In this direction, 


m = 4—>-H = 5.15T and therefore, due to the huge value 
of vf, the electric field that allow a remarkable effect is 
of the order of 10 6 V/m. Thus, figure [2] presents both 
contributions, S per (T , B, E) and S osc (T, B , E), as well as 
the magnetic entropy change AS(T, AB , E), considering 
E = 5 x 10 6 V/ m. Note the magnetic entropy change 
with an applied electric field AS(T, AB, E) is quite dif¬ 
ferent from the case without electric field A S(T, AB , 0). 
The reason to choose 104.5 K will be clear further in the 
text. We emphasize this temperature is quite higher in 
comparison to the case of a non-relativistic 3D diamag¬ 
netic material (c.a. 1 K) 7 . 

The influence of the electric field (E = 5 x 10 6 V/m), on 
the magnetic entropy change can also be seen as a func¬ 
tion of temperature, as presented in figure [3j The con¬ 
tributions S OS c(T, B, E) and S per {T 1 B 1 E) to the mag¬ 
netic entropy change AS(T, AB, E) are shown in fig¬ 
ures [31(a) and (b), for m = 1 and m = 2, respectively. 
Note S per (T, B, E) = 0 for m = 2, due to A m . Figure 
[3Kc) shows the influence of the electric field on the mag¬ 
netic entropy change, by comparing AS(T, AB, E) and 
A S(T, AB , 0), for some values of m. 

The behavior of the temperature dependence of the 
magnetic entropy change is rule by 7(z) function (see 
equation 1181) . that peaks at z max = 1.6 and therefore 
(from equation fT51) : 


T 

max 


^V~F N 0 


v~f N n 3B 2 


( 22 ) 


The first term only depends on the magnetic field and the 
second is quadratic on the electric field. Thus, an electric 
field decreases the temperature in which the magnetic 
entropy change is maximum (see figure [3](c)). At this 
point, it is possible to understand the temperature chosen 



FIG. 2. Magnetic entropy change with applied electric 
field AS(T, AB, E) (equation 1211) and its contributions 
S OS c{T, B, E) and S per {T, B, E). For the sake of comparison, 
also included is the case without applied field A S(T, AB, 0). 
These entropies are presented as a function of m, inversely 
proportional to the magnetic field (see equation 1121 . 
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to figure |H i.e., T = T max for m = 1 and E = 5 x 10 6 
V/m. 

On the other hand, an electric field changes the max¬ 
imum value of the magnetic entropy change (see figure 
|3](c)). To quantify the changing on the magnetic entropy 
change due to an applied electric field, let us consider an 
enhancement factor: 

= AS(T max ,AB,E) 

V AS(T max ,AB,0) 1 j 

&COS (' T max , B , E) + S s th(T max , B , E) 
S C os(T max , B , 0) 

For odd values of m: A m = cos(?u. 7 r) = — 1 and therefore 

V ~ 1 + n 2 L x ^/N 0 TT/3 - ^/3 2 (24) 

« 1 + 18.40/3 - 0.75/3 2 

Considering E = 5 x 10 6 V/m and B = 20.6 T (to = 
1), thus /3 = 0.2427. These values lead to 77 = 5.42, 
and means that an electric field increases the magnetic 
entropy change, since 77 > 1 (see figure [31(c)). On the 
other hand, for even values of to: A m = 0 and cos(to7t) = 


1. Thus 

V ~ 1 “ ( 25 ) 

Considering E = 5 x 10 6 V/m and B = 10.3 T (to = 2), 
thus /? = 0.4854. Note these values lead to 77 = 0.82 and 
therefore an electric field decreases the magnetic entropy 
change, since 77 < 1 (see figure E](c)). 

It is worth to note the magnetic entropy change for 
these results are not, as usual, for AB : 0 —> B, since the 
model does not allow B = 0. From equation [5J B m i n = 
E/vf and therefore B m i n = 5 T. 

Summarizing, we have described the influence of an 
applied electric field on the magnetocaloric properties of 
a 2D massless Dirac system: a graphene. For the case 
without electric field, this system presents an oscillating 
magnetic entropy change, due to a cosine term depend¬ 
ing on to, inversely proportional to the magnetic field. 
An electric field promotes an extra entropy term, with 
a periodic pattern also on to (this term has a maximum 
value for odd values of to and vanishes for even values of 
to). These two contributions give the magnetic entropy 
change, in which the maximum value either increases or 
decreases due to an applied electric field, depending on 
the value of to. Finally, the temperature in which the 
maximum magnetic entropy change occurs decreases due 
to an applied electric field. 
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(a) Magnetic entropy change A5(T, AB, E) 
and its contributions S 0 sc(T , E, E) and 
S P er(T , 5, E), for m = 1. 



T (K) 

(b) Magnetic entropy change AS(T , AE, E) 
and its contributions Sosc(T , .E) and 

Sper(T , E, E), for m = 2. 



(c) Comparison of A5(T, AE, E) and 
AS(T, A E, 0), for some values of m. It 
represents the influence of the electric field on 
the magnetocaloric effect of graphenes. 


FIG. 3. Temperature dependence of the magnetic entropy 
change of graphenes under a longitudinal applied electric field. 












